On the Absence of 0(a) Errors in Staggered-Quark Discretizations 
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We demonstrate that the O(a) taste mixing exhibited in standard textbook presentations of 
staggered quarks is an artifact of the particular definition of the flavor fields in those presentations, 
and has nothing to do with the underlying precision of staggered-quark actions, despite continuing 
comments to the contrary in the current literature. To illustrate this point we introduce a new 
coordinate-space definition of the flavor fields that suppresses the O(a) term by two additional 
powers of a. In fact there are no errors at all from this mechanism. The only source of taste mixing 
comes from the exchange of highly-virtual gluons and enters in 0(a 2 ). We review the idiosyncrasies 
of Symanzik improvement for naive/staggerd-quark actions, and show how these results follow from 
that program. 
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Improved versions of the staggered-quark action are 
proving highly effective for precise numerical simulations 
of full QCD pHa]- These formalisms have the unusual 
property that a single quark field creates several differ- 
ent flavors or "tastes" of quark, each with the same mass. 
This multiplication of quark degrees of freedom is un- 
physical but can be remedied to the extent that the differ- 
ent tastes decouple from each other • Unfortunately, 
despite almost thirty years of study, there remains con- 
fusion about the origins and size of taste mixing in these 
theories. It is critically important to resolve this issue 
given the central role played by taste mixing in the error 
analysis of staggered-quark results and given the theoret- 
ical issues raised by the presence of taste mixing 0-Q ■ 

Much of this confusion is caused by what have become 
standard textbook treatments of the flavor structure of 
staggered quarks [l(J. These treatments, which ignore 
interactions, re-express the free staggered-quark action 
in terms of a set of interpolating fields, with one field for 
each taste of quark. There are two standard choices for 
this set of fields: the coordinate-space flavor basis, and 
the momentum-space flavor basis. When rewritten in 
terms of the canonical coordinate-space flavor fields, the 
staggered-quark action has an 0(a) taste- mixing term, 
of the schematic form 



(1) 



where our notation is from The momentum-space 
basis, however, has no taste-mixing at all in the nonin- 
teracting case. Notwithstanding this discrepancy, sev- 
eral authors continue to assume that taste-mixing inter- 
actions occur in Ola). 

In fact, as we discuss in this note, the O(a) term in the 
coordinate-space analysis is highly misleading. It arises 
only because the coordinate-space flavor fields break one 
of the symmetries of the underlying staggered-quark the- 
ory: translation invariance on the original lattice. The 
0{a) term is needed to cancel errors introduced by the 
flavor fields. It has nothing to do with the underlying 
theory. Generally it is a bad idea to break the under- 



lying symmetries when forming interpolating fields, and 
there is no need to do so here. The momentum-space fla- 
vor basis, for example, does not violate translation sym- 
metry. Indeed the coordinate-space basis is never used 
in practical calculations; all such calculations are much 
more easily framed in terms of the momentum-space ba- 
sis . 

To underscore the artificial origins of the canonical 
0(a) term we will construct here a new coordinate- 
space flavor basis for which the taste-mixing term is 
0(a 3 ) rather than 0(a). Since the underlying theory 
is the same in both cases, neither the O(a) term nor 
the 0(a 3 ) term can be relevant to anything other than 
representations of the free theory in terms of the corre- 
sponding, flawed operators — in other words, neither is 
relevant to anything important. The momentum-space 
analysis gets the right answer: taste-mixing is impossible 
in the noninteracting case. As we will show, it is forbid- 
den by the symmetries of the underlying staggered-quark 
theories. 

Contrary to what is sometimes asserted, there is only 
one source for taste mixing and that is the exchange of 
highly virtual gluons between quarks [f| [l3l - [l5j | . Obvi- 
ously this mechanism cannot be analyzed in the noninter- 
acting theory, and therefore any analysis of the noninter- 
acting staggered-quark theory has nothing to say about 
real taste mixing. Here we will review what is known 
about real taste mixing. In particular we will review the 
arguments, some now quite old, for why these effects and 
all other finite-a errors enter first only in 0(a 2 ). We will 
also discuss how these errors show up in the spectrum of 
the staggered-quark Dirac operator, which is important 
for understanding vacuum polarization and for under- 
standing instanton effects. 

To address the 0(a) term, we only need examine 
free quarks since that term appears in the canonical 
coordinate-space analysis of the free theory. Rather than 
use the staggered-quark discretization directly, we will 
work with the formally equivalent but intuitively much 
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simpler "naive" discretization of the quark action 



with 



i>{x + aP)-^ X -_a^ + - {x)m ^ x) (2) 
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The gamma matrices and other conventions used here are 
described in Q. 

This action is invariant under a "doubling transforma- 
tion" 

ip(x) B ( (x)ip(x) $(x) -> ^(x)Bj(as) (3) 

where 

B c (x) = 7? (-l) Wa 

« n (T5^p) Cp exp(ix • Ctt/«): ( 4 ) 
p 

and C is a four-vector of Os and Is (i.e., £ M S Z 2 ) @. The 
"conjugate" £ to £ is (see Appendix A in (f|): 



Here we use the convenient notation 
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where n is a four-component vector with € Z 2 , to 
label the sixteen independent gamma matrices. 

This symmetry means that any low-momentum 
mode il>(x) is exactly equivalent to fifteen other modes, 
B((x)ip(x) for C 7^ 0. These modes are obviously dif- 
ferent from the original mode and from each other since 
Bq{x)iI){x) has very large momentum p w C 71 "/* 2 if ^{ x ) 
is a low- momentum mode (see Eq. (U)). Consequently 
the single naive-quark field actually creates sixteen dif- 
ferent flavors or "tastes" of quark. The field manages 
this trick by packing each of the different quark tastes 
into a different corner of its Brillouin zone, with taste £ 
corresponding to momenta p near (^ir/a. 

Sixteen is fifteen too many quarks, but the extra de- 
grees of freedom are easily removed if the different tastes 
do not mix (see @ for a review). The central issue, there- 
fore, is how much mixing there is between different tastes. 
The canonical analyses address this issue by re-expressing 
the action in terms of interpolating fields for each taste. 
In the canonical coordinate-space analysis, the lattice is 
divided into hypercubes that have two sites per side, and 
the interpolating fields are defined by weighted averages 
of the naive-quark field over the hypercube. The field 
corresponding to taste £ is [l2| 
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S B ( (x b )iP(xb) 
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where xb, with xb^/ci mod 2 = 0, identifies the hyper- 
cube, and 5 is a smearing operator defined by 



(8) 



S^tpix) 



ip(x) + ip(x + afi) 



(9) 



The smearing operator is critical to this definition. 
Its role is most easily understood in momentum space, 
where the operator becomes a simple function of mo- 
mentum 11611: 



S^p) = e^°/ 2 cos(p M a/2) 



— > 



1 for p^ — > 

ipfj.a/2 for p^ -)• it /a + p. 



(10) 
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It is used in the definition of \&«) to suppress p^ ~ ir/a 
while leaving p^ « unchanged. To see why this is 
important, we Fourier transform Eq. ([7]) to obtain 

* (C) (P)= X S(p + £r/o)7^(p+(C + 0T/o) (12) 



where we have used the Fourier transform of the doubling 
transformation Eq. ([3]), 



i/j(p) -» 7^ i>{p+ (ir /a). 



(13) 



Momentum p in \E^)(p) is restricted to — n /2a < p^ < 
■n /2a since the field is defined on the blocked lattice, 
which has lattice spacing 2a; this is also why we need 
the sum over £. The terms in ^^\p) with £ ^ cor- 
respond to tastes other than £, but they are suppressed 
by the smearing function which is 0(ap) unless £ = 0. 
Consequently, 

* (C) (p) -> l^^(p + C^/a) + 0(paijj) for \p^\ < ir/a. 

The smearing operator, working together with B^, serves 
to isolate momenta in the vicinity of Qir/a — that is, it 
isolates taste £, as is needed for the interpolating field. 

It is obvious from this analysis that the isolation of 
individual tastes is only approximate; VP^) is contami- 
nated by other tastes, £ + £, because the smearing op- 
erator suppresses these tastes by factors of only O(ap). 
Given that ty^, by its definition, contains other tastes, it 
seems likely that its field equation will exhibit 0(a) taste 
mixing terms. 

We can verify that this is the case by examining a 
solution ip{p) of the naive-quark Dirac equation, 



E 



i sin(»„a) \ , . , . 
^^7^ + m ip(p) = 0. (15) 



Note that this equation can be rewritten 

^(-1)C» is ™M lfi + m ^(p + C7r/a) = 0, (16) 
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or, equivalently, 

(£ + rn) 7? ^ + CTr/a) = 0, (17) 

which is a restatement of the doubling symmetry. 

Consider the blocked field with £ = 0, for example. 
It is easy to show that the Dirac equation, Eq. (TTSj) , for 
ijj(p) implies a blocked-lattice Dirac equation for ^^(p) 
with a small residue: 



i sin(2p M a 



^ 2a 



^7A. + m] *<°>(p) = A (18) 



where the residue A is (obviously) 

/ i sm(2 Ptl a) _ ^sin(p M a) 



x^ S (P + fr/a)ip(p + £7r/a). (19) 

Terms with £ ^ cause taste mixing, and when £^ = 
1 these are suppressed only to the extent that S(p + 
f7r/a) is small. In the continuum limit, |p M | <C 7r/a, the 
largest contributions to A therefore arise when only one 
component of £ is nonzero, £ = /t, and therefore, using 
Eqs. CEU) and QU, we find that 

A « ^(2i^)5(p + fi-K/a)ip(p + /i7r/a) (20) 



•^(-ira^75* (/i) (p). 

fi=0 



(21) 



As expected, A contains O(a) taste mixing. 

This 0(a) taste mixing is the conventional textbook 
result [l7j |. Our analysis shows that the suppression fac- 
tor of (ap) comes directly from the smearing operator. 
It is easy to design a new smearing operator for which 
the suppression is much stronger. For example, we can 
replace S M by 



S^tp(x) = 

ip(x — afi) + 3tp(x) + 3-i{j(x + ap.) + tp(x + 2ap) 



where now 



_ iP(to/2 3 cos(p M a/2) + cos(3 PA1 a/2) 



1 for p^ -> 

-i(p^a) 3 /8 for p M ->• n/a + p. 



(22) 

(23) 
(24) 



With this new smearing, the taste-mixing part of the 
residue A in the equation for the blocked field (Eq. (TT5)) ) 
is suppressed by order (pa) 3 rather than (ap): 

3 

A -> A = £ (-1)" (a^/4) 75 ^) (p) + ■ • • , (25) 

fj,=0 



where we have shown only taste-mixing terms. 

Obviously one can suppress taste mixing to even higher 
orders in the lattice spacing by using increasingly non- 
local smearings. Indeed one can suppress taste mixing 
completely, in the noninteracting theory, by simply re- 
defining 

¥<\p) = lx ij(p + C7r/a)H0(\p^\ < n/2a). (26) 



In coordinate space this would involve smearing over the 
entire lattice. 

This last definition is, of course, the standard 
momentum-space flavor basis. Using this basis it is triv- 
ially obvious that taste mixing is impossible in the free 
theory because it would violate momentum conservation. 
Momentum conservation is required by the translation 
symmetry of the underlying theory. Taste-mixing terms 
appear when we use the coordinate-space basis because 
that basis is defined using a hypercubic blocking of the 
lattice that breaks translation invariance and therefore 
violates momentum conservation, as is evident from the 
£ sum in Eq. (|12p . The momentum-space flavor basis, on 
the other hand, allows us to analyze separate tastes with- 
out violating momentum conservation and so is much 
more useful for studying finite-a errors in the theory. 

The impossibility of taste mixing in the free theory fol- 
lows immediately from the formal Symanzik analysis of 
finite-a corrections to the full interacting naive-quark ac- 
tion. Any finite-a correction to the original action must 
be a local operator that preserves all of the symmetries 
of the original theory, including gauge symmetry, trans- 
lation invariance, doubling symmetry, parity, and so on. 
The combination of the doubling and translation sym- 
metries implies that quark bilinears, in both the inter- 
acting and free theories, must have flavor-spin structure 
7n ® 1 [HI- For naive quarks these bilinears have the 
form 



where 



ipln <8> IV 1 i>(x) r y n ip(x ± Sx n ) 



6x!± 



(27) 



(28) 



guarantees the doubling symmetry. Consequently, these 
operators must be taste singlets in the noninteracting 
theory and cannot mix tastes; taste mixing is not allowed 
in the free theory! 

Such quark bilinears can and do mix different tastes 
when they are coupled to other fields, like the gluon field. 
This happens, for example, when the gluon field in a 
quark-gluon vertex carries off momentum q « £7r/a from 
the quark line. Such a momentum transfer leaves an on- 
shell quark still on shell, but with a different taste. The 
factor (— \)^' x / a in the gluon field alters the taste content 
of the quark bilinear: in general, an operator 



1p(x)^ n 1p(x ± Sx n ) X (-1) 



Q-x/a 



(29) 
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where the phase comes from a gluon (or other) field cou- 
pled to the quark bilinear, behaves in effect as though 
it had spin-taste structure j n -( ® £c 0- Consequently 
even a simple operator, like ip(x)j f j,U fl (x)ip(x+afl), which 
has nominal spin-taste structure 7^ x 1, actually has a 
very complicated spin-taste structure since the gluon field 
can carry off momenta with one or more components of 
order ir/a. 

The gluon in a taste-mixing interaction of this sort is 
highly virtual and so must be reabsorbed almost imme- 
diately by the same or another quark. Consequently the 
effects of such a taste-exchange interaction are indistin- 
guishable, to first approximation, from a taste-exchange 
four-quark operator. Since four-quark operators are at 
least dimension six, this means taste exchange can arise 
first in 0(a 2 ). Again this result follows immediately 
from the symmetry restrictions on quark bilinears in a 
naive/ staggered- quark theory: the only way a quark bi- 
linear in a naive-quark action can mix tastes is by cou- 
pling to a highly virtual gluon or gluons. This is the only 
mechanism for taste mixing allowed by the symmetries 
of naive/staggered-quark theories. 

The doubling symmetry obviously complicates the in- 
terpretation of correction terms in the quark action when 
interactions are included. Indeed it is more complicated 
still since the possibility that an on-shell quark can emit 
a highly virtual gluon and yet remain on shell means 
that additional factors of a<9 M acting on the gluon field 
in a quark bilinear do not necessarily suppress the op- 
erator by additional factors of a as would normally be 
the case [La ]. Despite these complications, the analysis 
of low-order Symanzik corrections to the naive-quark ac- 
tion is straightforward. In 0(a), the only operators that 
might enter are the standard 

mpl <g> lD 2 ip and mpgF -a® lip, (30) 

but each of these is suppressed by an additional factor 
of am, where m is the quark mass, because of the chiral 
symmetry of the naive-quark action when m — (see Q , 
for example). Consequently there are no 0{a) corrections 
of any sort; 0(a) errors are impossible for gauge-invariant 
quantities. 

The absence of 0(a) errors for physical quantities com- 
puted with staggered quarks has been known for a long 
time [l8j]. This should also be true of off-shell quanti- 
ties that are gauge invariant, such as eigenvalues of the 
naive/staggered-quark Dirac operator, since the suppres- 
sion of 0(a) errors does not rely upon redundancy. Re- 
dundant operators, which vanish by the equations of mo- 
tion, have no effect on physical quantities but can shift 



off-shell amplitudes. While some of the 0(a) operators 
are redundant, none of them contributes in 0(a) because 
they all violate the action's chiral symmetry and so are 
separately suppressed by an additional am. In the spe- 
cific case of the Dirac-operator eigenvalues, furthermore, 
it is hard to see how the eigenvalues could have 0{a) er- 
rors without causing similar errors in vacuum polariza- 
tion contributions to physical quantities, since the quark 
determinant in the path integral is the product of these 
eigenvalues. There is no way that such errors could be 
cancelled by valence-quark effects since these depend dif- 
ferently on the number of quark flavors rif. Consequently 
the absence of 0(a) errors in physical quantities requires 
that such errors are absent in the eigenvalues as well, as 
we expect. Recent simulation results for the eigenvalues 
also indicate that errors are likely 0{a 2 ) and smaller [20| . 

Strictly speaking, the Symanzik analysis is justified 
only for a straight simulation with naive quarks, includ- 
ing all sixteen tastes (or staggered quarks, with four 
tastes). After extensive study, all evidence thus far 
strongly suggests that the 0(a 2 ) taste- mixing in the 
sixteen-flavor theory translates into 0(a 2 ) errors in the 
"rooted" theory , where the quark determinant is re- 
placed by its l/16 th root in order to reduce the effective 
number of tastes to one. In particular no mechanism has 
been identified that could change an 0(a 2 ) error in the 
unrooted theory into an 0{a) error in the rooted theory. 
There is certainly no way that the 0{a) term from the 
canonical coordinate-space flavor basis (Eq. ([T])) could 
reappear since it violates momentum conservation and 
so was impossible to begin with; and there are no other 
unsuppressed 0{a) operators. 

The leading taste-mixing effects, therefore, are from 
the exchange of gluons with momenta of order Qn/a, 
where C,^ £ Z2. Luckily these, and all higher-order 
taste-mixing operators are highly perturbative, and so 
can be rigorously analyzed and systematically removed. 
We have already removed both the leading and next-to- 
leading order taste- mixing in going from naive/staggered 
quarks to ASQTAD quarks, and, more recently, from 
ASQTAD to HISQ quarks @|. The net effect has been to 
reduce taste-mixing effects by almost an order of magni- 
tude at current lattice spacings, as is immediately evident 
from nonperturbative calculations of the mass differences 
between pions of different taste [2l|. Taste mixing is 
clearly under control. 
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